A Two Potential Formula and its Application to Proton-Proton Scattering 
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for the on-shell partial wave scattering amplitude. This formula embodies a large number of possible 
applications, including long range Coulomb forces as well as short distance singular potentials. As 
an example illustrating the use of the formula we analyze the determination of the strong proton- 
proton scattering s-wave phase shift from the experimentally determined Coulomb phase when the 
one pion exchange and two pion exchange chiral potentials are taken into account and analyze the 
relevant scales of the problem. 
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I. INTRODUCTION 

The two-potential formalism, developed in the fifties 
by Watson [l[ and Gell-Mann and Goldberger 0, re- 
lates the scattering due to the sum of two different po- 
tentials and has a widespread use in scattering theory. 
The usual example is the treatment of Coulomb distor- 
tion for strongly interacting particles. The problem is to 
determine the total scattering amplitude T from a poten- 
tial constructed as the sum of two potentials V — Vs + Vl 
in terms of the scattering amplitude T$ due only to the 
potential Vs- The result can be found in a straightfor- 
ward manner using the Lippmann-Schwinger equation 



T = V + VGqT, 



(1) 



with V the potential operator and Gq = (E — i?o) _1 the 
resolvent of the free Hamiltonian. The outgoing bound- 
ary condition corresponds to E — > E+i0 + . The T-matrix 
can then be expressed as 1 



T = T s + (1 + T S G ) f L (G T S + 1) 



(2) 



where Ts is the short distance T-matrix, and Tt, the long 
distance one distorted by short range effects 



V s + VsG Ts , 
V L + V L G S T L , 



(3) 
(4) 



with Gs — Go + GqVsGs the full propagator for Vs- 



While the result above solves the problem, it does not 
explicitly relate the on-shell scattering amplitudes of the 
full and short distance potentials. The reason is that the 
Lippmann-Schwinger equation involves the off-shell be- 
haviour of the potentials, which allows to treat nonlocal 
potentials quite straightforwardly, but precisely because 
of this it is hard to profit specifically from the simplifying 
features which arise in the interesting and quite frequent 
case of local potentials arising e.g. in a particle exchange 
picture. For the local case, a coordinate space formu- 
lation of the scattering problem is more convenient (an 
effective field theory example is provided by Ref. Q). 

In this paper we derive a two-potential formula, which 
relates the phase shifts (i.e. the on-shcll scattering ma- 
trix) of the full and short range potentials V and Vs, 
and which is based on two assumptions: (i) the poten- 
tials are local and (ii) the short range potential domi- 
nates at short distances. Our result will be amenable 
to rather detailed analytical study, hence enlarging the 
class of situations one may cover. The connection to mo- 
mentum space renormalization with counterterms of the 
Lippmann-Schwinger equation is also analyzed. This is 
particularly enlightening in the case of singular poten- 
tials and their renormalization, a subject of recent inter- 
est (see e.g. ||). As an illustrative application we show 
how our two-potential formula may be used to deduce 
the strong proton-proton phase shifts from the experi- 
mentally measured ones when long distance corrections 
coming from one and two pion exchange contributions 
are taken into account. 
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1 We have chosen this particular formulation of the two-potential 
trick in order to have a more consistent notation along the paper. 
Normally it is written in terms of the distorted short range T- 
matrix and the undistorted long range T-matrix, but in cither 
case the resulting T-matrix is the same. 



II. TWO POTENTIAL FORMULA 

We consider the non-relativistic scattering of two par- 
ticles by a spherically symmetric potential V which can 
be decomposed into the following two pieces 

V(r) = V L (r) + V s {r), (5) 



2 



where Vl and Vs respectively represent the long and 
short distance components of the interaction. We will 
assume that the short range potential Vs dominates at 
short distances r — r c , i.e. 



V s (r B )>V L (r c ), 



(6) 



for r c small enough. The system can be described by 
solving the reduced Schrodinger equation (for simplicity, 
we only consider here the s-wave case) 

- u'l + lix [V L (r) + V s (r)} u fc (r) = k 2 u k (r) , (7) 

with u k (r) the reduced wave function, fj, the reduced mass 
of the two body system and k the center-of-mass momen- 
tum. The phase shift can then be obtained by matching 
the reduced wave function u k to the usual asymptotic 
boundary condition for r — > oo 



Uk (f) — ► cot 5 sin kr + cos kr . 



(8) 



where we have assumed that the long range potential Vl 
decays faster than 1/r 2 at large distances, so phase shifts 
are well defined. We also consider the corresponding scat- 
tering problem for which only the short range potential 
Vs is present. In such a case, the reduced Schrodinger 
equations reads 



u b k +2»V s (r)u b k (r) = k A u b k (r) 



(9) 



with u k the short reduced wave function. The phase shift 
can be extracted from the asymptotic behaviour of u k 



u ki r ) ~ > c °t $ sm kr + cos kr , 



(10) 



for r — ► oo. 

The problem is to relate the full phase shift 8(k) with 
the short phase shift 6 s (k). For that purpose, we will 
assume that at short enough distances r — r c the full and 
short reduced wave function u k and are approximately 
equal, u k (r c ) ~ u k (r c ). The previous approximation can 
be restated in terms of the logarithmic derivatives of the 
reduced wave functions 2 



u 'kiTc) 
u k (r c ) Uk( r c) 



4'( r c 



(11) 



This expression will hold true when the condition ex- 
pressed in Eq. © is fulfilled. We now make use of the 
superposition principle to represent the full solution as 
the following linear combination 



u k {r) = cot S J fc (r) - Y k (r) , 



(12) 



where Jk and Y k are solutions of Eq. ([7]), subjected to 
the asymptotic conditions 



Mr) 
Yk(r) 



sin kr , 
— cos kr . 



(13) 
(14) 



2 We are not necessarily assuming a regular solution of the 
Schrodinger equation, i.e. Vg can contain zero-range pieces. 



Analogously, we write the short range solution as 

u s k (r) = cot 5 s J£(r)-Y k s (r), (15) 



with and Y k s solutions of Eq. ([9]), such that 



4 S W 



Y k S (r) 



sin kr , 
— cos kr . 



(16) 
(17) 



for r — > oo. By matching the logarithmic derivatives we 
arrive at our final expression 



cot S(k) 



A(k,r c ) cot 6 s (k) - B{k,r c ) 



C(k,r c ) cot S s {k) -V{k,r c 
where A, B, C and T> are defined as 

A(k,r c ) = W(jS,Y k )\ r=rc 
B(k,r c ) = W(Y k s ,Y k )\ 



C(k,r c ) 
V(k,r c ) 



W{JtJ k )\ r=r 
W(Y k s ,J k )\ 



(18) 

(19) 
(20) 
(21) 
(22) 



with W(f,g)\ r = rc = f(r c )g'(r c ) - f'(r c )g(r c ) the Wron- 
skian between different wave functions evaluated at the 
cut-off radius r — r c . The bilinear structure is reminis- 
cent of the Moebius transformation invariance discussed 
at length in Ref. [J] in the context of the renormalization 
group analysis with boundary conditions. It should be 
noted that although the matching of log-derivatives in 
order to obtain long range correlations is not new, one 
nice example being the Landau-Smorodinsky derivation 
of the effective range expansion or the treatment of 
hadronic atoms in Ref. [6j, its use in combination with 
the superposition principle in order to derive direct re- 
lations between phase shifts is less common, and it has 
only been partially exploited in some effective field theory 
r-space computations @, d, Hi • 

In passing we also note that Eq. (TT8"|) cannot be de- 
rived from the Lippmann-Schwinger equation. The rea- 
son is that the two-potential formula depends on the 
explicit formulation of the following: (i) the superposi- 
tion principle via Eqs. (TT2"|) and (Ti"5"|) . and (ii) the short 
distance boundary condition for the Schrodinger equa- 
tion, Eq. (jlip . These two conditions are included in 
the Lippmann-Schwinger equation, but implicitly, in a 
way that they cannot be directly handled, impeding 
the derivation of the previous formula (but allowing the 
derivations of formulas relating the full off-shell scatter- 
ing amplitudes, like the two-potential trick, Eq. ((2|)). 

The two potential formula also holds in certain cases 
for non-local potentials. The necessary condition for its 
application is that the non-local potential does not in- 
volve derivatives of order higher than two, e.g. potentials 
of the type 



V 



NL 



{V 2 ,/s(0}, 



(23) 



where {,} represents the anti-commutator. In such 
a case the short distance boundary condition for the 
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Schrodinger equation can be expressed as the log- 
derivative of the wave function, Eq. (fTTj) . For non-local 
potentials involving higher derivatives, the two-potential 
formula can still be applied when the cut-off radius r c is 
larger than the range at which the non- localities appear. 



III. COULOMB SCATTERING 

The case where the long range potential is of Coulomb 
type requires a special treatment as the usual asymptotic 
behaviour described in Eq. ([8]) does not apply. For defi- 
niteness, we analyze here the Coulomb repulsion between 
two unit charge particles in the s-wave. The full system 
is now described by the following equation 



,C 



2p 



V s (r) + - 
r 



u c k {r) = fcXf(r), (24) 



2„,CV 



where we have added the c superscript for labelling the 
Coulomb solution, and a represents the fine structure 
constant. The correct asymptotics for u k is given by 



u£(r) -> cot5° F (i 1 ,p) + G ( V ,p), 



(25) 



with 5° the Coulomb modified phase shift, and F (r],p) 
and Go(r],p) the usual s-wave Coulomb wave functions 
(see for example [lfj|V which depend on the parameters 
r\ = l/(ka B ) and p = kr; a B = l/(pa) is the Bohr radius 
of the two particle system. Fo and Go are solutions of the 
reduced Schrodinger equation for the Coulomb potential 
Vc(r) = a/r, with the asymptotic behaviour 

F (r],p) — > sin (kr - r/log2kr + a ) , (26) 
G {il,p) -> cos(fcr-77log2fcr + cr ), (27) 

where <tq is a phase shift defined as Co = argT(l + irj). 
As in the previous case, we can use the superposition 
principle to rewrite the full solution 

u£(r) = cot5 c F^{r)-G c k {r), (28) 

with Fk(r) and Gfc(r) solutions of Eq. (|24|) subjected to 
the asymptotic boundary conditions 

FC(r) -> F (v,p), (29) 
G^(r) -> -G {v,p). (30) 



The short range system is described by Eq 

,s 



and the 



short range wave function u k is again parametrized by 



(31) 



Eq. dT5J). 

After matching logarithmic derivatives we find 

c A(k,r c ) cot 6 s (k) -B{k 7 r c ) 

C ° [ ' C{k,r c ) cot 5 s (k) -V(k,r c ) 

where A, B, C and T> are now defined as 

A(k,r c ) = W(J^G c k )\ r=rc 

B(k,r c ) = W(Y k s \G^)\ r=rc 

C(k,r c ) = W{Jl,F£)\ r=rc 

V(k,r c ) = W(Y k s ,FC)\ 



(32) 
(33) 
(34) 
(35) 



in complete analogy with Eq. (jT8|) . Previous relationships 
between Coulomb and short distance scatte ring can be 
found for some specific cases in Refs. [TT1. fl2l. fl3| . 



A. Contact Short Range Potential 

A simple application of the previous formula corre- 
sponds to a situation where the short range potential 
is zero for distances greater than the cut-off radius r c 



Vs(r) = for r > r c 



(36) 



while, for distances shorter than r c , the potential is very 
strong. The previous potential corresponds to a 5-type 
contact interaction regularized at the length scale r c . In 
such a case, the F k (r) and G k (r) wave functions are 
equal to their asymptotic behaviour for r > r c , and by 
taking into account their behaviour at small radii 



G C k (r) 



k C{ri) 
1 



0(r 3 



-0(r 2 ) 



r 

a B 

2r / 2r 

log — 

a B V a B 



2 1E - 1 



(37) 
h( V ) 
(38) 



with je the Euler-Mascheroni constant, and C(rj) and 
h(rj) defined as 



c 2 (v) 



2TT7] 



(39) 



Hv) = V 2 Y] f 2 , 2 > - log?? - IE , (40) 
* — ' n(n + r> 

n—l v ' 

the relationship given by Eq. (|3"Tj) can be evaluated ex- 
plicitly, yielding 



kcotb s (k) = C 2 {rf)k cot 8 c (k) + 2 



a B 



2 

a B 



a B 



l 0g _iL_2 7£ +0(r c ), (41) 



where terms linear in the cut-off radius and higher pow- 
ers of r c have been ignored. As can be seen the previ- 
ous expression is logarithmic divergent with respect to 
r c , but can be regularized if we take into account the 
corresponding expression for k = 0, which is similar to 
the well known relationship between stron g an d Coulomb 
scattering length from Blatt and Jackson jl4l. [l5j 



1 

as 



1 

ac 



2 

a B 



a B 

log— -2 1E 



+ 0{r c ). (42) 



The expression above diverges in exactly the same way as 
Eq. (|4Tj) . Subtracting the k = expression to Eq. (|4*Tj) . 
and taking the r c — > limit, we arrive at the following 
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expression 



k cot 5 s (k) + — 

a s 



C 2 (j])k cot S c (k) 



a B 



etc 



(43) 



The expected error of this formula can be estimated rein- 
troducing the cut-off r c and interpreting it as the ne- 
glected range Rs of the short range potential Vs, yielding 
a relative error of 0{r c /aB) = 0(Rs/cib)- 

The corresponding formula for attractive Coulomb in- 
teraction may be of interest for the treatment of pionic 
atoms, and can be obtained by taking r\ = — l/ikas) 
negative, and making the following substitution 



h(rj) -> Rc [i'(iT]) - log {-irf)} 



(44) 



with ip the digamma function. 

Finally, the corresponding formula for p-wave repul- 
sive Coulomb interaction can be worked out analogously 
to the s-wave case. The treatment of the divergences is 
nonetheless more involved: there is an additional loga- 
rithmic divergence proportional to fc 2 , due to the inter- 
play between the Coulomb potential and the centrifugal 
barrier. The outcome is that two subtractions are needed 
in order to have finite results in the r c — > limit. The 
final formula is rather simple to summarize 



fc 3 cot5f(fe) 



1 



1 



r hS k 2 = C*(rj) fc 3 cot <jf (jfe) 



"1,5 

+ fc 3 (l 
1 



r l 2 )2r,h{r t ) 
1 



T,n,ck 2 



(45) 



where C 2 (i]) = (1 + r] 2 )C 2 (ri), oti t s and c*i,c are the p- 
wave short and Coulomb scattering volumes, and n,s 
and r\fi the p-wave short and Coulomb effective ranges. 
The previous formula has less predictive power than the 
corresponding one for s-waves as a consequence of the 
extra subtraction needed to regularize it. A possible ap- 
plication is nucleon-alpha scattering [l6[ . 



IV. APPLICATION TO PROTON-PROTON 
SCATTERING 

Now we apply the previous results for Coulomb scat- 
tering to the specific case of proton-proton (pp) scatter- 
ing in s-waves. We consider the strong pp interaction as 
the short range potential Vs, while the Coulomb repul- 
sion between the protons plays the role of the long range 
potential Vl- 



A. Pionless theory 

We first consider the simplifying case in which the pion 
exchange interactions between the protons are neglected 



and the pp potential consists on contact interactions only, 
i.e. the pionless theory, characterized by a short distance 
boundary condition. In such a case the two potential 
formula given by Eq. (|43[) applies. The previously men- 
tioned relationship can be better understood by noticing 
the relationship with the strong and Coulomb effective 
range expansions [13], i-e. 



k cot S b = - 



1 



rs k 2 



a s 2 

OO 

$>„, S fc 2n , (46) 



71 = 2 



k cot S c C 2 (rj) + 2 



Hv) 



1 1 



a c 2 

OO 



rc k 2 



(47) 



n=2 



meaning that, with the exception of the scattering length, 
which explicitly depends on the regularization scale r c , 
see Eq. (|42"T) , the strong and Coulomb effective range pa- 
rameters for pp scattering are equal in the present ap- 
proximation 



rs = r c and v n s = v n , c for n > 2 , 



(48) 



where rs(c) is the effective range, V2,s(c) the shape pa- 
rameter, etc. If we compare the previous results with the 
parameters obtained with the Nijmegen II potential [l8l ]. 
we observe a small discrepancy 3 



r 5 = 2.84fm , r c = 2.76fm, 



(49) 



giving a 3% relative difference between the strong and 
Coulomb parameters. According to the error estimation 
of the previous section, we should expect a relative error 
of Rs/dB, with Rs the range of the strong pp interac- 
tion, given by one pion exchange, Rs = R^o = l/m n o, 
where m T o the neutral pion mass, yielding the result 
R^o/olb = M p a/2m^o ~ 2.5% (M p is the proton mass) 
in agreement with the previous discrepancy 4 . 

The corresponding results for the strong pp phase 
shifts, obtained from the Coulomb pp ones for the Ni- 
jmegen II potential [lH by means of Eq. (j43|) . are shown 
in Fig.[TJ The agreement between the strong pp Nijmegen 
II phase shifts and the expected ones computed from the 



Instead of the Nijmegen II values, it is also possible to use the 
well-established experimental value for the Coulomb pp effective 
range rc = 2.794(14) fm, and the model dependent strong pp 
one rs = 2.84(4) fm (see Ref. [3), although the conclusions do 
not change appreciably. 

The contributions to the strong-Coulomb effective range dif- 
ference from vacuum polarization [20ll . or from the modified 
Coulomb potential of Ref. [2l| , are expected to be much smaller 
than the strong (pionic) corrections. The magnetic moment in- 
teraction [22l | does not contribute for s-wave proton-proton scat- 
tering. 
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FIG. 1: (Upper panel) Strong pp phase shifts computed from 
the Coulomb pp phase shifts (Nijmegen II potential [l8|) by 
using the zero range strong-Coulomb correlation of Eq. (|43[) 
and its corresponding extensions when including the OPE 
and chiral TPE (N 2 LO) potentials with a cut-off radius of 
r c = 0.1 fm. (Lower panel) Difference between the Nijmegen 
phase shifts and those obtained with the strong- Coulomb cor- 
relations. 



contact theory correlation (143)) is quite good, never ex- 
ceeding a 1.5° difference, as expected from the relative 
error estimation. 



B. Comparison with other approaches 

The above result may be relevant to the effective field 
theory (EFT) formulation of low-energy pp scattering 
done by Kong and Ravndal [HI, [Hj] based on the power 
divergence subtraction (PDS) regularization scheme of 
Refs. [25|, l26(- The admitted intricacy of the momen- 
tum space formalism in those works contrasts with the 
much shorter and transparent discussion of the coordi- 
nate space renormalization presented above. In partic- 
ular, Eq. ([4"3")l implies that there are no Coulomb cor- 
rections to the effective range once the cut-off is re- 



order calculation of Kong and Ravndal 24], but dis- 
agrees with the next-to-next-to-leading order results of 
Ref. [27J • It is also implied in our treatment that the 
pionless treatment of pp scattering can be made almost 
scale independent if, apart from the usual counterterms 
Co + C*2(p 2 + p' 2 ) + 0(p 4 ), a counterterm contribution 
proportional to a is included in the computations, i.e. 
aD e 2. This observation is closely related with the re- 
sults of Ref. [28[, where the necessity of a strong and 
Coulomb version of Cq was discussed. The previous D e 2 
counterterm fixes the difference between the strong and 
Coulomb scattering length, so the price to pay in order 
to remove the log scale dependence in Eq. ([12")) is the im- 
possibility to relate the two, as both scattering lengths 
become input parameters. This seems to be in contra- 
diction with Kong and Ravndal [HI, [24| , who argue that 
the C2 counterterm stabilizes the scattering length (see 
also related discussions in Refs. [27], |29|, |30| ) . This coun- 
terterm is analogous to the D2 counterterm needed 
to renormalize Weinberg power counting at leading or- 
der [3l], H3] ■ They are both due to the similar behaviour 
of the Coulomb and Yukawa potential at short distances. 

Of course, these conclusions are based on our coor- 
dinate space analysis with cut-off regularization. Di- 
mensional regularization with PDS yields different re- 
sults as Coulomb corrections to the effective range 
appear at next-to-next-to-leading order. These correc- 
tions depend on the off-shell behaviour of the 0(p 4 ) coun- 
terterms, which in the cut-off approach seems to be un- 
der control as long as non-localities and off-shell ambi- 
guities happen below r c . A more pessimistic view is pre- 
sented by Gegelia in Ref. [2§ |. where it is argued from 
the renormalization group behaviour of the counterterms 
in PDS that it needs to be a strong and Coulomb ver- 
sion of each counterterm (to absorb the log-divergences) , 
meaning that in the end it is impossible to relate strong 
and Coulomb observables. On the contrary, the renor- 
malization group analysis with cut-off regularization of 
Birse and Barford [33j seems to support the idea that 
the Coulomb log-divergences can be absorbed in just one 
counterterm 5 . The observations of Gegelia [28| can be 
considered as an extension to any scattering observable 
of the results of Refs. [H, [Hj] about the difficulty of 
obtaining model independent strong scattering lengths 
from Coulomb ones due to short range ambiguities of the 
wave function. The previous seem to be in contradiction 
with usual requirement of short distance independence 
of physical results in effective field theory. In fact, as 
was shown in Ref. (36j, further constraints on the short 
range ambiguities not considered in [3J, |35j can notice- 



5 The results of Ref. [33l | does not exclude the existence of 
Coulomb corrections to all counterterms, neither do our results 
if supplemented by additional subtractions. It is just that they 
are not needed in order to have scale independent results. 
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ably reduce the model dependence of strong parameters, 
in a better agreement with EFT expectations. Finally, 
we should also stress that we are only trying to separate 
strong from Coulomb corrections in non-relativistic quan- 
tum mechanics. A complete formulation on the separa- 
tion of strong and electromagnetic effects is only possible 
in the context of quantum field theory, see Ref. [33] for a 
modern discussion on the subject. 

We should nonetheless remember that cut-off regular- 
ization is a physical regularization, in the sense that the 
cut-off r c can be interpreted as a physical scale. From this 
point of view, the meaning of the log-divergence in the 
relationship between the strong and Coulomb scattering 
length is straightforward: it represents the expected error 
of the strong scattering length in the pionless approxima- 
tion, which scales as \og(Rs/aB) (instead of Rs/ob, as 
in the other parameters), yielding a very large, about 
~ 350%, expected relative error (to be compared with 
the one for the effective range ~ 2.5%). One can also ar- 
gue that the extra counterterm D e i is not needed, as the 
inclusion of the higher order components of the potential 
will reduce the scale dependence. 

C. Strong-Coulomb Correlations with Chiral TPE 
Potentials 

The two potential formula makes it possible to ob- 
tain the (experimentally inaccessible) strong pp phase 
shifts from the (experimentally accessible) Coulomb ones. 
While a complete analysis should of course include vac- 
uum polarization [2(| , modifications to the Coulomb po- 
tential [2l| and even exchange [H, HiJ H(| , the in- 
teresting issue is that one can obtain model independent 
strong phase shifts, provided that we employ a model 
independent strong pp potential Vs and model indepen- 
dent Coulomb phase shifts. Here, we will do so with 
an eye put on the relevant scales in the problem, an as- 
pect which our two-potential treatment can address in a 
rather clean way. 

For the previous purpose we use the potentials of chiral 
perturbation theory [4j| as the short distance potential 
V s of Eq. |[3"T]). These potentials include TPE effects and 
can be expressed as a expansion in powers of Q 

V s (r) = V® (r) + V<® (r) + V™ (r) + 0(Q 4 ) , (50) 

where Q represents either the pion mass or the momen- 
tum of the protons. We also use the Nijmegen PWA 
which is a model-independent extraction of the pp s-wave 
Coulomb 6 phase shifts from a large proton-proton scat- 
tering database. With that information, we can obtain 



It is important to notice that the pp phase shifts in the Nijmegen 
PWA are not Coulomb phase shifts, but electromagnetic phase 
shifts. By that it is meant that the pp phase shifts are defined 
with respect to the asymptotic solutions of the full electromag- 
netic potential used in their analysis, which consists on improved 



the strong pp phase shift and its error from the PWA 

6° WA (k) ± A^ WA -> (V x , r e ) -> 5 s (k) ± A5 S , (51) 

and analyze the resulting cut-off dependence, which is 
an important issue, as for large coordinate space cut-offs 
the higher order pieces of the chiral potential are not 
resolved. 

It should be noted here that a complete model inde- 
pendent separation between strong and electromagnetic 
contributions is not always possible, specially if short dis- 
tance electromagnetic effects arc included. One exam- 
ple is the inclusion of nucleon form factor corrections to 
the magnetic moment interaction in the proton-proton 
PWA of Ref. [43| . Another example is proton finite size 
corrections to the Coulomb potential. The formalism 
presented here clearly separates between what we define 
as the strong and electromagnetic potential. That does 
not necessarily mean that exact model independence has 
been achieved, specially if corrections like the ones men- 
tioned above are added, or that strong and electromag- 
netic effects have been actually separated, specially as 
electromagnetic corrections to the proton mass or to the 
coupling constants have not been included. 

The specific procedure we will apply is analogous to 
the one followed in the pionless case, i.e., we do not 
directly use the strong- Coulomb two-potential formula, 
Eq. (|31[) . but rather perform a subtraction of the equiva- 
lent two-potential formula for the scattering lengths, and 
then check for cut-off independence of the results. This 
choice also allows for a better comparison between the 
pionless correlation given in Eq. (|43|) and the correspond- 
ing improvements when the strong physics are included 
explicitly. 

In the present calculation we are only going to consider 
the chiral potentials up to the Q 3 order, i.e. next-to-next- 
to-leading-order (N 2 LO). At this order the finite range 
piece of the chiral potential consists of one-pion exchange 
and chiral two-pion exchange. An interesting feature of 
the chiral two-pion exchange potentials is that they are 
highly singular, diverging as ~ 1/r 6 at N 2 LO 7 . In har- 
mony with previous findings @, [H, Q , this divergence will 
become rather unimportant: the two potential formula 
shows a smooth cut-off dependence for singular chiral po- 
tentials 8 . In any case, we stress that our main concern is 
to analyze the minimal short distance cut-off r c for which 



Coulomb, vacuum polarization and the magnetic moment inter- 
action (see Ref. 14211 for details). As our current analysis is not 
intended to be complete, we will ignore most of these details and 
simply assume that the long range potential is the usual Coulomb 
potential, and that the full electromagnetic phase shifts roughly 
coincide with the Coulomb ones, <5p WA — <5pwA- 

7 The most singular (non-contact) piece of Vi will behave as 
l/r 3 +" in coordinate space and as \q\" f (\q] / m-^) in momentum 
space, being q the momentum exchanged between the nucleons 
and / a non-polynomial function. 

8 In fact, the singular chiral two-pion exchange potentials yield 
smoother results than the OPE potential. In the current regu- 
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FIG. 2: Strong pp s-wave phase shifts at £a a b = 50 MeV and 
-Blab = 200 MeV as a function of the cut-off r c , c omp uted from 
the Coulomb pp phase shifts (Nijmegen PWA [421 ]') by using 
the chiral potential truncated at different orders in the chiral 
expansion (pionless means no potential, OPE is the leading 
order potential and TPE the full N 2 LO chiral potential). The 
strong and Coulomb phase shifts have been related by making 
use of Eq. (|31[1 with one subtraction at k — 0, in order to have 
finite results when the cut-off is removed. 



higher order effects can be distinguished from lower order 
ones, rather than the specific cut-off dependence of the 
results. 

The results for E lah = 50 MeV and E Xah = 200 MeV 
can be seen in Fig.O By TPE we refer to the N 2 LO chiral 
potential, and for compactness we skip the NLO results, 
which lie between the LO (OPE) and the N 2 LO results. 
The bands represent the error coming from the original 
pp Coulomb phase shift in the PWA of Ref. [Hj]. For the 



larization scheme, OPE shows a mild log-divergence at distances 
of 10 — 3 fm. This divergence can be eliminated by using a dif- 
ferent, and more complex, subtraction prescription, but for the 
purposes of the present discussion it is not particularly important 
what happens at such small scales. 



strong and Coulomb scattering lengths, which are needed 
for the subtractions, we take the values corresponding to 
the Nijmegen II potential [Hj], i.e. ac = — 7.81 fm and 
a 5 = — 17.25fm. The actual expressions for the chiral 
pp potential are taken from Ref. [4l[. Only the long 
range piece of the chiral potentials is considered, and 
the corresponding counterterms are ignored: they are 
equivalent to a boundary condition for the Schrodinger 
equation [H, [45[ , and are therefore already implicitly in- 
cluded in the two-potential formula. We take = 1.29, 
m^o = 134.98 MeV and f n = 92.4 MeV, which accord- 
ing to the definitions of Ref. [42j gives an / 2 = 0.0755 
for the scaling mass m s = m K ± = 139.57 MeV. The 
previous chiral pp potential explicitly depend on three 
chiral couplings, ci, C3 and C4, which appear at 0(Q 3 ) in 
the expansion of the potential, and which relate nucleon- 
nucleon and nucleon-pion scattering. We take for these 
chiral couplings the values obtained in Ref. [46| from an- 
alyzing the pp data alone, i.e. c\ — —0.76(7) GeV -1 , 
c 3 = -4.78(11) GeV -1 , and c 4 = 3.92(52) GeV -1 . As 
can be seen, for a cut-off above r c — 1.2 fm, one can- 
not distinguish, within uncertainties, between lower and 
higher order computations, i.e. it does not matter 
whether pions are included in Vs or not. Actually cut- 
offs below r c = 0.8 fm are needed in order to fully distin- 
guish the chiral two-pion exchange contributions within 
the accuracy of the phase shifts. This result is not en- 
tirely surprising as could be anticipated from considering 
the two-pion exchange related Compton wavelength scale 
^2tt = l/2m,r ~ 0.7 fm. A more complete analysis should 
include vacuum polarization, 2irj and 77 exchange ef- 
fects, which will affect the the precise values of the strong 
phase shifts but will hardly change the observation on 
the relevant scales. The same remarks also apply to the 
error analysis, which should include the error in the sub- 
tracted strong and Coulomb scattering lengths and the 
theoretical uncertainties in the chiral potential itself, like, 
for example, the error in the determination of the chiral 
couplings. 

We can also compare the extracted strong effective 
ranges for the different cases considered. In this case, 
it is used the Coulomb pp phase shift from the Nijmegen 
II potential as input for the two-potential formula and 
the resulting strong phase shifts are shown in Fig. [TJ For 
the regularization scale r c = 0.1 fm, we obtain 

rs.contact = 2.78 fm, (52) 
rs,OPE = 2.63 fm, (53) 
rs,TPE = 2.87 fm, (54) 

to be compared with the Nijmegen II result, r$ = 2.84 fm. 
The pionless value differs from the one given in Sec- 
tion I IV Al due to finite cut-off effects, while OPE sur- 
prisingly contributes in the wrong direction. The TPE 
result reproduces the Nijmegen II one within a 1% accu- 
racy level, and agrees within error estimations with the 
extraction of Ref. [l!|, rg = 2.84(4) fm, where the error 
accounts for different sources of model-dependence. 
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Finally, we note that despite the TPE potential be- 
comes highly singular at short distances, diverging as 
~ 1/r 6 , nothing dramatic happens, making the limit 
r c — > innocuous precisely when the TPE effects become 
visible, i.e. for r c < 0.8fm 9 . This particular feature is a 
specific merit of our two potential formula which provides 
a clean separation between scales and implements in an 
extended distorted wave fashion the renormalization pro- 
gram carried out in previous works (see e.g. [3, 



V. CONCLUSIONS 

The two potential formalism provides a framework 
where forces of different origin and ranges may be dis- 
entangled rather explicitly. We have proposed a coor- 
dinate space formulation which re-states the result in a 
rather transparent way and fully exploits the boundary 
value character as well as the superposition principle of 
the scattering problem. Our result allows for a detailed 
investigation of the relevant scales built into the problem. 
This is particularly enlightening in the case of singular 
potentials and their renormalization, a subject of recent 
interest. We have exemplified our approach by discussing 
its consequences for the proton-proton system, where the 
electromagnetic and strong forces contribute to the scat- 
tering process, as a method to extract the strong phase 
shifts in a model independent fashion. We have only 
discussed s— waves (with the exception of Eq. (f4"5|) ) and 
single channel scattering. The extension to higher par- 
tial waves, as well as coupled channels, is straightforward 
but cumbersome, see Appendix|A] Such an extended for- 
malism might allow to discuss further interesting appli- 
cations of the present ideas to similar problems where a 
scale separation of different forces would be necessary. 
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APPENDIX A: EXTENSION TO HIGHER 
PARTIAL WAVES AND COUPLED CHANNELS 

1. Higher Partial Waves 

The extension of our two potential formula to higher 
partial waves is straightforward. The full two-body 
system is described by the corresponding reduced 
Schrodinger equation for the Z-wave 



<i + 



2^(V L (r)+V s (r)) + fc^ 



Uk,i(r) 



= k 2 u kil (r), (Al) 

where Uk,i(r) is the /-wave reduced wave function, \x the 
reduced mass and k the center-of-mass momentum. The 
asymptotics of Uk,i for r — > oo is given by 



Uk,i(r) -> cot (5/ ji(kr) - yi(kr) , 



(A2) 



where ji and yi{kr) are the reduced spherical Bessel func- 
tions, defined as ji(x) = x ji{x) and yi(x) = xyi(x). We 
only consider here the case of a long range potential Vl 
decaying faster than 1/r 2 at large distances. By mak- 
ing use of the superposition principle we rewrite the full 
solution as 



Uk.i (r) = cot Si Jk,i (r) - Y k ,i (r) , 



(A3) 



where Jk,i and Yk,i are solutions of Eq. (|A1[) subjected 
to the asymptotic boundary conditions J kt i{r) — > ji{kr) 
and Yfc,;(r) — > yi(kr) for r — » oo. The corresponding 
scattering problem for which only the short range poten- 
tial Vs is present is described by the reduced Schrodinger 
equation 



' k.I 



2/iVs(r) 



1(1 + 1) 



= k 2 u s k ,(r) 



(A4) 



with u? , the short l-w&ve reduced wave function. The 
short distance phase shift is obtained from the asymp- 



totic behaviour of u k , 



u li( r ) -> cot Sfjt(kr) - yi(kr) , 



(A5) 



for r — > oo. We rewrite itf l as 



9 Actually, despite Eq. (0 being fulfilled to the extreme, the TPE 
correction is small because the typical minimal wavelength is still 
not comparable to the range where the TPE correction takes 



ul l (r)=cot6 l Jl l (r)-Y k s l (r) : 



(A6) 



with J kl and Y kl solutions of Eq. (|A4[) obeying the 
asymptotic boundary conditions JuAr) — > ji(kr) and 
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As usual we match the logarithmic derivatives of 
Uk,i( r ) an d Uki( r ) a t the cut-off radius r — r c , yielding 



for r — > oo, where j and y are diagonal matrices given by 



Ai(k,r c ) cot<Sf(fc) -g t (fc,r e ) 
COt6lW C^k, r c ) cot af (Jfe) -2Mfc,»0 

where Si, C/ and I?; are defined as 

Mk,r c ) = w(j^,y fc) ,)| r=ro , 
s ; (/c ; r c ) = ^(Y£,Y*,i)| r=rc > 

A(/c,r c ) = W(^,J M )I =r , 



(A7) 



(A8) 
(A9) 
(A10) 
(All) 



in analogy with the s-wave case. In principle the use of 
the previous formula is straightforward as long as a finite 
cut-off is employed in the computation. On the contrary, 
if one tries to remove the cut-off, some divergences may 
appear, mostly related to the centrifugal barrier. There- 
fore a detailed analytical or numerical study of the diver- 
gences will be necessary in order to obtain a finite result 
in the r c — > limit. 



2. Coupled Channels 

The extension to the coupled channel case is direct to 
obtain if an adequate notation is used. We will consider 
the general case of N coupled channels. They can be 
described by the following Schrodinger equation, which 
in compact notation reads 



2/z(V L (r)+V s (r)) + T 



Ufc(r) 



= fc 2 u fe (r-), (A12) 

where the wave function u& is now an N x N matrix, 
each column representing a linearly independent solu- 
tion. The long and short range potentials and V5 
are also N x N matrices (the non-diagonal terms relating 
the different channels), and L 2 is the angular momentum 
matrix, which is diagonal and given by 



diag(Ji(li + l)Mh 



1) 



,MZiv + l)),(A13) 



being l\, I2, In the orbital angular momentum of each 
channel. In principle there are 2N linearly independent 
solutions (two per channel), but regularity conditions at 
the origin reduce this number to N. This is why the 
wave function can be represented by an N x N matrix. We 
have also added the simplifying assumption that there are 
no inelasticities, meaning that in Eq. ([AT2j) the source 
of the coupling is cither tensor forces or dipole-dipole 
interactions. 

The asymptotic behaviour of the wave function matrix 
Ufc is given by the following expression 



j(fcr) = diag(j (l (kr),j h (At), . 
y(fcr) = dia,g(y h (kr),yi 2 (kr), 



,jl N (kr)), (A15) 
,y, N (fcr)),(A16) 



with ji and yi the reduced spherical Bessel functions as 
defined in the previous section. The matrix M(fc) is the 
analogous of cot 5 for coupled channels and is related to 
the S-matrix by M(fc) = i (S(fc)+1) /(S(fc)-l) with 1 the 
identity matrix. It is a symmetric matrix and contains 
N(N + l)/2 independent scattering parameters or phase 
shifts. By making use of the superposition principle, we 
can rewrite the wave function matrix as 



u*(r)=J fc (r)M(fc)-Y fc (r), 



(A17) 



where J& and Y& are solutions of Eq. (j A12|) which asymp- 
totically behave as Jfe(r) — ► j(fer) and Yfc(r) — * y(kr). 

The corresponding Schrodinger equation for the short 
wave function is 



2/xV s (r) + 



ug(r) = k 2 u s k (r) , (A18) 



where, in analogy to the full case, the short wave function 
matrix can be written as 

uf(r)=jf(r)M s (fc)-Yf(r), (A19) 

with J|? and Yf solutions of Eq. (|A18[) subjected to 
the asymptotic boundary condition Jf (r) — > j(Ar) and 
Yf(r) — > y(kr) for r — > 00. 

For obtaining the corresponding two potential formula 
one needs to match the logarithmic derivatives of the 
wave functions, which for the coupled channel case means 

UfcVcXu^rc))- 1 = uf (r c )(uf (r c )) _1 . (A20) 

Using the following Wronskian relation 

ul(r c )u k '(r c ) = uf (r c )u fe (r c ), (A21) 

where the T superscript denotes the transpose, the 
boundary condition given by Eq. (|A20|) can be rewrit- 
ten as 

uj(r c )uf(r c )=u2"(rc)uf(r c ), (A22) 

an expression which does not involve the inverse of the 
wave functions. If we rewrite and in terms of 
Eq. (|A17[) and (|A19|) . we arrive at our final expression 

M{k) - (A(k,r c )M s (k)-B(k,r c )) x 

(C(fc,r c )M s (fc)-P(fc,r c )) _1 (A23) 



ujt(r) -> j(fcr)M(fc) -y(Ar) 



(A14) 



10 The reason why we write Jj.(r)M(fc) instead of M(fc) (r) 
in Eq. (I A 1 71 ) is because if is a solution of the Schrodinger 
equation ||A12| | and A a constant N X N matrix, then A is 
also a solution (but this is not the case for Au/j). 
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with A, B, C and V defined as where the Wronskian is given by W(F, G)| r = 

F'(r c )G(r c ) - F(r c )G'(r c ). 



A(k,r c ) = 


-T^(Yj,j£)| r=rc , 


(A24) 


B{k,r c ) = 


-W(Y^Y s k )\ r=rc , 


(A25) 


C{k,r c ) = 


-W(3l,3 s k )\ r=rc , 


(A26) 


V(k,r c ) = 




(A27) 
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